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Variations  in  the  Southern  Hemisphere  (SH)  winter  of  2007,  2008  and  2009  had  important  con¬ 
sequences  on  polar  mesospheric  clouds  (PMCs)  observed  in  the  corresponding  Northern  summers. 
Specifically,  the  stratospheric  SH  winter  of  2007  was  observed  to  be  warmer  than  in  2008  and  2009. 
Using  the  high  altitude  analysis  from  the  Navy  Operational  Global  Atmospheric  Prediction  System- 
Advanced  Level  Physics  High  Altitude  (NOGAPS-ALPHA)  forecast/assimilation  system  we  show  that  this 
warmth  was  linked  to  similar  temperature  increases  in  the  high  latitude  summer  mesosphere.  These 
temperature  changes  led  to  a  dramatic  reduction  in  PMC  occurrence  (factor  of  5-6)  recorded  by  the 
SHIMMER  instrument  at  sub-arctic  latitudes  and  a  factor  of  2  decrease  in  total  ice  water  content  in 
PMCs  seen  by  the  SOFIE  instrument  on  the  NASA  AIM  satellite.  Microphysical  modeling  confirms  the 
overall  effect  of  these  temperature  changes  on  PMCs  at  high  latitudes;  however,  a  detailed  comparison 
of  the  cloud  occurrence  with  the  SHIMMER  data  for  all  three  years  shows  that  the  clouds  are  associated 
with  a  surprisingly  wide  range  (130-165  I<)  of  temperatures. 
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1.  Introduction 

It  is  now  generally  understood  that  the  climatology  of  the 
middle  atmosphere  differs  greatly  between  Northern  and  South¬ 
ern  Hemispheres.  Due  to  greater  land  sea  contrasts  and  greater 
topography,  the  North  experiences  more  wintertime  planetary 
wave  activity,  which  leads  to  greater  wave  induced  heating  in  the 
NH  winter  stratosphere,  a  weaker  and  more  easily  perturbed  NH 
polar  vortex,  and  reduced  breaking  of  gravity  waves,  which 
contribute  to  the  warm  winter  stratopause  (Yulaeva  et  al.,  1994; 
Garcia  et  al.,  1992;  Hitchman  et  al.,  1989).  In  recent  years,  interest 
has  shifted  to  what  these  differences  might  mean  for  the  summer 
seasons  in  each  hemisphere.  Early  work  by  Alexander  and 
Rosenlof  (1996)  showed  that  the  summer  stratosphere  is  warmer 
in  the  SH  relative  to  the  NH  due  to  greater  gravity  wave  induced 
forcing  in  the  southern  summer.  Siskind  et  al.  (2003)  showed  that 
these  stratospheric  hemispheric  asymmetries  had  mesospheric 
counterparts  whereby  there  would  be  weaker  gravity  wave  drag 
in  the  Southern  upper  mesosphere.  This  was  suggested  as  a 
possible  cause  of  the  reduced  occurrence  of  polar  mesospheric 
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clouds  (PMCs)  in  the  southern  summer  since  weaker  gravity  wave 
drag  would  imply  a  warmer  summer  mesopause  (Siskind  et  al., 
2005a). 

Most  recently  there  has  been  great  interest  in  so-called 
teleconnections  between  wintertime  conditions  in  one  hemi¬ 
sphere  and  the  corresponding  summer,  which  occurs  simulta¬ 
neously  in  the  other.  For  middle  atmospheric  conditions,  this 
interest  was  provoked  by  the  unusual  Southern  winter  of  2002. 
During  this  winter,  several  minor  stratospheric  warmings  were 
observed  (Siskind  et  al.,  2005b),  culminating  in  an  unprecedented 
major  warming  in  September  (e.g.  Allen  et  al.,  2006  and 
references  therein).  At  the  same  time,  the  MACWAVE  rocket 
campaign  (Goldberg  et  al.,  2004)  reported  unusual  conditions  in 
the  NH  summer  mesosphere.  These  included  a  warmer  meso¬ 
pause  and  reduced  occurrence  of  polar  mesospheric  summer 
Echoes.  Model  simulations  by  Becker  et  al.  (2004)  and  Becker 
and  Fritts  (2006)  led  to  the  suggestion  that  the  unusually  strong 
planetary  wave  activity  in  the  Southern  Hemisphere  initiated  a 
chain  of  events,  which  ultimately  led  to  perturbations  to  the 
gravity  wave  drag,  which  produces  the  cold  summer  mesopause. 

Subsequent  work  by  Karlsson  et  al.  (2007,  2009a  and  2009b) 
have  shown  that  these  long  range  teleconnections  between  the 
winter  in  one  hemisphere  and  the  summer  in  the  other  occur 
quite  often.  This  effect  can  be  documented  for  both  interannual 
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variability  as  well  as  within  a  single  season.  In  general,  the  greater 
planetary  wave  activity  in  the  NH  winter,  and  its  large  inter¬ 
annual  variability,  has  a  major  effect  on  the  year-to-year  varia¬ 
bility  of  SH  PMCs.  However,  even  though  the  effect  is  weaker, 
a  connection  between  the  weaker  waves  in  the  SH  winter  and  NH 
PMCs  could  also  be  detected.  The  Karlsson  et  al.  (2007)  analysis 
relied  on  a  meteorological  analysis,  which  extended  only  up  to 
1  mb,  thus  excluding  the  mesosphere. 

Here  we  use  an  analysis,  which  extends  up  to  90  km  to  look  at 
interannual  variability  in  NH  PMC  from  2007-2009.  Section  2 
presents  3  NH  seasons  of  mesospheric  cloud  data  from  two 
instruments,  the  Solar  Occultation  for  Ice  Experiment  (SOFIE) 
and  the  Spatial  Heterodyne  Spectrometer  for  Mesospheric  Radi¬ 
cals  (SHIMMER).  SOFIE  is  on  board  the  NASA  Aeronomy  of  Ice  in 
the  Mesosphere  (AIM)  satellite  (Russell  et  al.,  2009),  launched 
in  April  2007  and  as  of  this  writing,  still  making  measurements 
of  the  high  latitude  region  where  PMCs  form.  SHIMMER  was 
the  primary  payload  of  STPSat-1,  launched  in  March  2007  for  a 
mission,  which  lasted  30  months  (Englert  et  al.,  2010).  Both 
SHIMMER  and  SOFIE  observed  mesospheric  clouds  during  the 
2007,  2008  and  2009  Northern  summer  seasons,  but,  as  we  will 
discuss,  using  quite  different  observational  techniques  and  from 
different  orbits.  Section  3  describes  results  from  the  synoptic 
analysis  provided  by  the  NOGAPS-ALPHA  forecast/assimilation 
system.  In  Section  4,  we  describe  the  use  of  NOGAPS-ALPHA 
temperature,  water  vapor  and  winds  to  drive  the  Community 
Aerosol  and  Radiation  Model  for  Atmospheres  (CARMA)  micro¬ 
physical  model.  We  calculate  ice  water  content  and  compare  that 
with  SOFIE  observations.  Finally,  Section  5  concludes  and  sum¬ 
marizes  the  major  findings  of  this  work. 


2.  Observations  of  mesospheric  clouds:  SOFIE  and  SHIMMER 

2.1.  SOFIE 

SOFIE  uses  the  solar  occultation  technique  to  retrieve  vertical 
profiles  of  PMC  extinction  at  eleven  wavelengths  ranging  from  the 
near  UV  to  the  mid  IR  (Gordley  et  al.,  2009).  A  salient  feature  of 
SOFIE  is  that  its  observations  are  confined  to  a  single  latitude 
band.  Since  AIM  is  in  a  polar  sun-synchronous  orbit,  this  latitude 
varies  in  a  narrow  range  throughout  the  season,  increasing  from 
65°  in  late  June  to  near  70°  in  August  (cf.  figure  3  of  Gordley  et  al., 
2009).  Also  SOFIE  sampling  remains  at  a  single  local  time,  near 
2300  LST  for  the  Northern  summer  PMC  season. 

Hervig  et  al.  (2009a)  gave  an  overview  of  data  from  the 
northern  summer  of  2007.  They  noted  that  after  July  1st,  there 
were  periodic  episodes  of  PMC  decreases  coincident  with  periodic 
warmings  of  the  summer  mesopause  region  of  up  to  6  K.  Fig.  1 
compares  the  frequency  of  occurrence  for  the  first  three  Northern 
seasons.  The  SOFIE  data  version  is  1.022.  Fig.  1  shows  that  SOFIE  is 
so  sensitive  that  it  often  sees  clouds  100%  of  the  time.  This  is 
particularly  true  for  2008  and  2009.  For  2007,  there  are  periods 
after  July  1st  where  the  cloud  occurrence  frequency  dips  to  as  low 
as  80%.  Hervig  et  al.  (2009a)  comment  on  these  reductions  in 
frequency  and  column  ice  mass  as  being  consistent  with  a 
warming  event  which  caused  PMCs  to  sublimate.  Fig.  1  shows 
that  these  mid-season  cloud  reductions  were  not  pronounced  in 
2008  or  2009. 

As  described  by  Hervig  et  al.  (2009b),  SOFIE  PMC  extinctions  at 
IR  wavelengths  are  directly  proportional  to  ice  volume  density, 
which  given  suitable  assumptions  for  ice  density,  can  yield  the  ice 
mass  density.  The  vertical  integral  of  the  ice  mass  density  yields 
the  column  ice  mass,  also  known  as  ice  water  content  (IWC). 
Stevens  et  al.  (2010)  and  Siskind  et  al.  (2007)  have  shown  that  the 
IWC  is  a  useful  quantity  for  comparisons  between  microphysical 


Fig.  1.  Frequency  of  PMC  occurrence  observed  by  SOFIE  for  the  NH  summers  of 
2007,  08  and  09.  The  horizontal  axis  covers  the  period  from  about  May  21st  to 
September  21st.  The  latitude  of  the  SOFIE  occultation  roughly  tracks  the  sun  and  is 
about  66°  at  solstice,  rising  to  near  80°  in  September  (e.g.  http://sofie.gats-inc. 
com/sofie/index.php). 


Days  Relative  to  Solstice 

Fig.  2.  PMC  ice  water  content  measured  by  SOFIE  for  the  NH  summers  of  2007, 
2008  and  2009.  The  x-axis  is  the  same  as  Fig.  1. 


models  and  satellite  measurements.  Fig.  2  shows  the  SOFIE  IWC 
for  the  3  years  in  question.  Similar  to  the  occurrence  frequencies, 
starting  around  July  1  (10  days  relative  to  solstice),  the  IWC  in 
2007  stands  out  as  being  lower  than  in  2008  or  2009.  This  relative 
reduction  persists  until  early  August  (around  +45  days  relative  to 
solstice). 

2.2.  SHIMMER 

The  SHIMMER  data  nicely  complement  the  SOFIE  data.  While, 
like  SOFIE,  SHIMMER  viewed  clouds  on  the  limb,  instead  of 
occultation,  it  viewed  the  scattered  sunlight  from  the  clouds.  As 
discussed  in  Stevens  et  al.  (2009),  solar  scattered  light  from  PMCs 
occasionally  enhances  the  Rayleigh  scattered  signal  observed  by 
SHIMMER  near  309  nm  and  between  about  80  and  85  km  altitude. 
To  identify  PMCs  we  fit  an  exponential  function  to  the  back¬ 
ground  away  from  altitudes  where  PMCs  could  contaminate  the 
fit.  The  imposed  threshold  for  cloud  detections  is  furthermore  a 
function  of  the  solar  scattering  angle  to  minimize  the  effects  of 
lighting  variations  on  our  detections.  As  a  result  of  the  STPSat-1 
orbit,  SHIMMER  records  data  over  a  wide  range  of  local  times, 
with  a  precession  rate  of  about  1  \2  h/day.  Over  the  course  of 
the  PMC  season,  SHIMMER  sampled  all  local  times  for  which  the 
clouds  were  illuminated. 
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Because  of  the  low  inclination  of  the  STPSat-1  orbit  (37°), 
SHIMMER  only  observed  up  to  58°N.  The  data  we  showed  comes 
from  latitudes  between  50  and  58 °N,  roughly  corresponding  to 
the  latitude  region  where  ground  based  observers  see  these 
clouds  as  noctilucent  clouds  (NLCs).  We  also  only  considered 
clouds,  which  occurred  at  altitudes  above  80  km.  This  largely 
eliminates  the  potential  problem  of  clouds  occurring  in  the 
foreground  or  background  of  the  SHIMMER  line  of  sight,  which 
would  be  then  erroneously  ascribed  to  temperatures  that  might 
be  too  high  (i.e.  at  too  low  altitude).  This  consideration  is 
important  for  the  comparison  of  SHIMMER  cloud  occurrences 
with  NOGAPS-ALPHA  temperatures  in  Section  3. 

Fig.  3  shows  the  number  of  PMCs  observed  daily  by  SHIMMER 
for  the  three  Northern  seasons  (2007,  2008  and  2009)  that  it 
collected  data  (this  is  Version  2  of  the  SHIMMER  data).  Since  the 
daily  sampling  rate  was  nearly  identical  for  each  of  the  three 
seasons,  these  histograms  are  nearly  proportional  to  occurrence 
frequency.  Also  shown  is  the  local  time  of  the  observation,  which 
for  a  given  date,  is  slightly  different  for  the  three  years.  Since 
diurnal  variations  in  PMC  frequency  are  known  to  be  significant, 
care  should  be  taken  to  compare  the  three  seasons  at  as  close  a 
local  time  and  date  as  possible.  With  this  small  caveat  in  mind,  it 
is  readily  apparent  that  after  solstice,  the  number  of  clouds 
observed  in  2007  is  much  less  than  in  2008  or  2009.  For  example, 
for  11-12  LST  (near  July  13th-15th  in  2007),  less  than  10  clouds 
were  seen  whereas  in  2008  (near  July  8th-10th)  over  40  were 
seen  and  about  10-20  detections  in  2009  (near  July  2nd-4th). 

By  dividing  the  cloud  detections  by  the  total  number  of 
SHIMMER  observations  along  the  orbital  track,  we  calculate  a 
frequency  of  occurrence.  Table  1  presents  these  values  along  with 
monthly  averaged  occurrence  frequencies  and  IWC  for  SOFIE  for 
the  three  years.  The  table  shows  that  in  July  2008  and  2009,  SOFIE 
detected  a  cloud  on  98  or  99%  of  its  orbits.  In  July  2007,  the 
averaged  occurrence  frequency,  while  still  high,  was  less  at  92%. 
The  SHIMMER  frequency  is  less  than  1%  in  2007,  and  increases 
by  about  a  factor  of  5  to  6  in  2008  and  2009.  The  much  smaller 
occurrence  frequency  and  the  relatively  greater  change  in  fre¬ 
quency  seen  by  SHIMMER  relative  to  SOFIE  is  probably  partially 
due  to  the  reduced  sensitivity  of  SHIMMER  to  dim  clouds 
compared  with  SOFIE.  But  it  is  undoubtedly  also  due  to  the  fact 
that  SHIMMER  views  clouds  at  the  equatorward  edge  of  their 
occurrence  where  they  are  particularly  sensitive  to  small  changes 
in  the  background  saturation  conditions.  The  monthly  averaged 
SOFIE  ice  water  content  shows  that  the  July  2007  reduction  in 
PMC  occurrence  corresponded  to  a  factor  of  2  reductions  in  ice 
water  content.  Other  data  taken  using  ground-based  lidar  at  69°N 
(e.g.  Fiedler  et  al.,  2011)  show  an  analogous  dip  in  both  frequency 
and  brightness  in  2007  relative  to  adjacent  years,  at  least  for  the 
brightest  clouds  (which  tend  to  track  the  total  IWC) 


3.  Meteorological  overview 

To  quantitatively  understand  the  interannual  variability  of  the 
summer  mesopause  we  use  the  NOGAPS-ALPHA  forecast/assim¬ 
ilation  system,  specifically  the  synoptic  analysis  first  discussed  by 
Hoppel  et  al.  (2008)  and  Eckermann  et  al.  (2009).  The  NOGAPS- 
ALPHA  analysis  uses  the  NRL  Atmospheric  Variational  Data 
Assimilation  System  (NAVDAS).  NAVDAS  is  a  3D-variational 
(3DVAR)  system  with  a  6  h  update  cycle  that  assimilates  both 
conventional  operational  meteorological  data  as  well  as  data  from 
the  NASA  AURA  Microwave  Limb  Sounder  (temperature,  ozone 
and  H20)  and  TIMED  Sounding  of  the  Atmosphere  with  Broad¬ 
band  Emission  Radiometry  (SABER)  temperature  data.  Hoppel 
et  al.  (2008)  presented  the  first  realistic  mesospheric  analysis  and 
showed  results  for  January  and  February  2006.  They  assimilated 
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Fig.  3.  Frequency  of  PMC  occurrence  observed  by  SHIMMER  for  the  NH  summers 
of  2007,  2008  and  2009.  The  solid  line  is  the  number  of  PMC  observations  on  a 
given  day,  referenced  to  the  left  axis.  The  black  circles  are  the  associated  daily 
averaged  local  times  of  the  cloud  detections  and  the  orange  dots  are  the  daily 
averaged  local  times  of  all  the  SHIMMER  images,  regardless  of  whether  there 
were  clouds  or  not  (both  local  times  are  referenced  to  the  right  axis).  The  dates  for 
each  year  are  given  by  the  upper  x-axis.  Note  that  2008  was  a  leap  year  so  the 
dates  differ  by  1  day  relative  to  2007  and  2009.  (For  interpretation  of  the 
references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version 
of  this  article.) 


Table  1 

PMC  occurrence  frequencies  (%)  and  ice  water  content  (pg  m-2). 


Instrument 

Month 

2007 

2008 

2009 

SOFIE  (%,  pg  nr2) 

June 

95,  38 

95,  31 

95,  37 

SOFIE 

July 

92,  36 

99,  69 

98,  61 

SOFIE 

August 

74,  23 

82,  30 

75,  22 

SHIMMER  (%  only) 

July 

0.6 

3.2 

3.6 

2016 
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MLS  data  up  to  p  =  0.01  hPa  (about  80  km  in  geometric  coordi¬ 
nates).  Eckermann  et  al.  (2009)  extended  this  capability  by 
assimilating  new  versions  of  the  MLS  (Version  2.2)  and  SABER 
(Version  1.07)  data  up  to  0.002  hPa  (about  88  km)  and  produced 
an  analysis  for  the  summer  of  2007  to  support  the  NASA  AIM 
mission.  The  forecast  model,  which  is  run  concurrently  with 
NAVDAS  extends  up  to  p  =  0.0005hPa  (about  92  km);  we  will 
show  results  up  to  90  km.  Eckermann  et  al.  (2009)  also  presented 
some  validation  of  the  planetary  waves  and  tides  that  are 
captured  by  the  analysis  as  well  as  the  discussion  of  the  tuning 
of  the  gravity  wave  drag  parameterization.  The  Eckermann  et  al. 
(2009)  configuration  was  subsequently  frozen  as  a  production 
configuration,  and  used  to  generate  6  h  analyses  beyond  July 
2007,  out  to  (currently)  early  2010.  It  therefore  covers  the  three 
NH  PMC  seasons  under  consideration,  including  the  entire  dura¬ 
tion  of  the  SHIMMER  mission.  Stevens  et  al.  (2010)  used  an 
average  of  the  June  2007  analysis  fields  to  initialize  the  CARMA 
microphysical  model  and  to  calculate  the  diurnal  variation  of  PMC 
ice  water  content. 

Fig.  4  summarizes  the  zonal  mean  temperature  variability  in 
these  three  seasons  by  comparing  2008  and  2009  with  2007  for 
the  months  of  June,  July  and  August.  It  shows  that  the  tempera¬ 
ture  differences  in  June  were  small,  except  that  June  2007 
was  slightly  colder,  than  June  2008  and  2009.  This  is  consistent 
with  the  variation  presented  in  figure  9  of  Stevens  et  al.  (2010).  It 
is  also  qualitatively  consistent  with  the  fact  that  the  cloud 
frequency  in  the  SOFIE  data  differed  little  amongst  the  three  June 


months.  However,  for  July  and  August,  the  comparison  of  2007 
with  2008  and  2009  reveals  interesting  global  anomaly  patterns. 
The  Southern  hemisphere  shows  a  vertically  alternating  pattern 
of  warm  and  cold  anomalies.  This  pattern  is  similar  to  that 
recently  presented  by  Azeem  et  al.  (2011,  submitted  to  Journal 
of  Geophysical  Research).  They  present  ground  based  tempera¬ 
tures  over  Antarctica  that  suggest  anomalous  conditions  in  2007, 
probably  due  to  enhanced  levels  of  planetary  wave  activity  in 
2007  relative  to  2008.  The  tri-level  pattern  of  alternating  warm 
and  cold  anomalies  also  resembles  that  seen  during  the  SH  minor 
stratospheric  warming  observed  in  August  2002  and  reported 
by  Coy  et  al.  (2005)  and  Siskind  et  al.  (2005b),  although  the 
magnitude  of  the  stratospheric  temperature  increase  here 
(  +  15  I<)  is  much  less  than  the  50  K  reported  for  the  2002  event. 
The  anomaly  pattern  is  also  offset  from  the  southern  pole.  This 
feature  was  also  previously  commented  on  by  Karlsson  et  al. 
(2007)  as  consistent  with  the  weaker  wave  activity  in  the  SH 
relative  to  the  NH  causing  the  variations  to  occur  on  the  vortex 
edge  rather  than  at  the  highest  latitudes,  as  it  does  in  the  NH. 

Most  significantly  for  the  discussion  here,  Fig.  4  shows  that  the 
2007  stratospheric  warm  anomaly  appears  to  extend  upward  and 
northward  to  the  summer  upper  mesosphere.  This  is  also  very 
similar  to  the  patterns  identified  by  Karlsson  et  al.  (2009a)  in  the 
Canadian  Middle  Atmosphere  Model;  thus  warm  winter  strato¬ 
spheric  anomalies  correlate  with  warm  summer  mesopause 
anomalies.  The  summertime  morphology  is  similar  to  Karlsson 
et  al.  (2009a)  in  that  the  warm  anomaly  tilts  upward  and 
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Fig.  4.  Comparison  of  monthly  mean  zonal  average  NOGAPS-ALPHA  temperature  deviations  for  3  Northern  summer  months,  in  all  cases  referenced  to  2007. 
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poleward  into  the  arctic  mesopause  region.  Thus  at  sub-arctic 
latitudes,  the  entire  upper  mesosphere  is  warm;  however,  as  one 
moves  to  the  pole,  the  warm  anomaly  rises  in  altitude,  closer 
to  the  mesopause.  Finally,  the  July  2007  anomalies  are  larger 
than  the  August  2007  ones.  This  is  also  qualitatively  consistent 
with  the  interannual  PMC  variation  recorded  by  SOFIE. 

Since  PMC  observations  are  typically  referenced  to  actual 
geometric  altitudes,  it  is  useful  to  convert  the  NOGAPS  tempera¬ 
tures  to  an  altitude  grid  (cf.  Stevens  et  al.,  2010).  Furthermore, 
because  NOGAPS  provides  a  complete  synoptic  analysis,  it  is 
straightforward  to  sample  the  analysis  at  the  times  and  locations 
of  the  SOFIE  and  SHIMMER  measurements.  Fig.  5  shows  this  for 
the  latitudes  and  local  times  of  the  SOFIE  observations  along  with 
the  actual  SOFIE  temperatures  for  July  and  August  of  the  three 
years  in  question.  This  comparison  serves  as  a  valuable  compar¬ 
ison  of  the  NOGAPS  temperatures  with  an  independent  dataset. 


While  one’s  eye  is  drawn  to  the  difference  for  altitudes  above 
85  km  between  NOGAPS  and  SOFIE,  given  uncertainties  in  the  C02 
and  O  densities  required  by  the  SOFIE  non-LTE  retrieval,  this 
difference  may  not  be  surprising.  In  addition,  there  may  be  a 
cold  bias  to  the  NOGAPS  temperatures  because,  as  discussed  by 
Schwartz  et  al.  (2008),  MLS  appears  to  be  up  to  10  K  too  cold  at  the 
highest  altitudes  in  polar  summer.  As  noted  by  Eckermann  et  al. 
(2009),  we  apply  a  globally  averaged  bias-correction  to  the  MLS 
data;  this  correction  equals  4.3  K  at  0.002  hPa  (the  top  data  level). 
This  leaves  a  remaining  cold  bias,  which  in  principle,  could  equal 
5-6  K;  still  much  smaller  than  the  almost  20 1<  difference  between 
SOFIE  and  NOGAPS  in  Fig.  5.  However,  for  purposes  of  validating 
the  structure  of  the  temperature  changes,  we  argue  that  it  is  more 
useful  to  highlight  the  important  similarities  between  the  two 
datasets.  This  is  emphasized  in  Fig.  6,  which  compares  the  year-to- 
year  relative  changes  in  NOGAPS-ALPHA  and  SOFIE. 
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Fig.  5.  Monthly  mean  zonal  mean  temperatures  from  SOFIE  compared  with  NOGAPS-ALPHA  temperatures  sampled  at  the  same  latitudes  as  SOFIE  (see  the  caption 
of  Fig.  1). 


Fig.  6.  Comparison  of  temperature  differences  for  2008  and  2009  relative  to  2007  for  NOGAPS-ALPHA  (solid  lines)  and  SOFIE  (dashed  lines).  The  dotted  lines  are  presented 
just  for  reference  to  zero  temperature  change. 
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Fig.  6  shows  that,  at  all  altitudes,  even  where  they  disagree, 
both  NOGAPS  and  SOFIE  show  the  same  year-to-year  variability 
and  both  show  that  above  80-82  km,  2007  was  warmer  than 

2008  and  2009.  Further,  both  NOGAPS  and  SOFIE  also  agree  that 
the  warm  anomaly  was  greater  in  July  2007  relative  to  2008  and 

2009  than  it  was  in  August  2008  and  2009.  At  90  km,  July  the 
temperature  enhancement  in  2007  is  in  the  range  of  9  to  11  K  in 
both  SOFIE  and  NOGAPS;  in  August,  both  datasets  indicate  a  3-6  K 
temperature  enhancement.  This  is  strong  validation  of  the 
NOGAPS-ALPFIA  temperature  variability  near  its  top  boundary 
and  also  strongly  validates  both  the  global  pattern  of  the  tele¬ 
connection  and  its  detailed  vertical  structure  in  the  high  latitude 
summer  mesopause  region  suggested  in  Fig.  4. 

Also  of  interest  in  Fig.  6  is  that  the  interannual  temperature 
change  is  quite  small  at  cloud  altitudes  (82-83  km),  typically 
between  0  and  1 .5  K.  Flowever,  it  grows  rapidly  with  altitude  and 
at  85  km  it  is  in  the  range  2-5  K.  As  we  will  see,  it  appears  that  it 
is  the  temperature  difference  at  85  km,  which  is  important  for  the 
greater  ice  abundance  in  2008  and  2009. 

At  the  latitudes  sampled  by  SHIMMER,  no  SOFIE  data  is 
available  so  we  just  show  mean  temperatures  from  NOGAPS  for 
July  of  the  three  years  in  question  (SHIMMER  did  not  detect  PMCs 
in  August).  Fig.  7  shows  that  the  overall  temperatures,  which 
correspond  to  the  SHIMMER  observations  are  much  warmer,  and 
the  2007  warm  anomaly  is  greater,  than  at  SOFIE  latitudes.  At 
90  km,  2007  temperatures  exceed  2008  and  2009  by  15-20  K. 
The  mean  temperatures  of  150-160  K  are  too  high  to  form  PMCs. 
This  is  consistent  with  the  much  lower  frequency  of  occurrence  at 
SHIMMER  latitudes  relative  to  the  SOFIE  latitudes.  While  SOFIE 
sees  PMC  almost  all  the  time,  PMCs  are  a  relatively  rare  event  in 
the  SHIMMER  database,  generally  associated  with  excursions 
below  the  mean  temperature. 

Fig.  8  illustrates  the  effect  of  selecting  only  those  temperatures 
in  NOGAPS-ALPHA  associated  with  SHIMMER  cloud  observations. 
It  presents  two  histograms  of  SHIMMER  temperatures  for  the 
three  years.  This  is  similar  to  Fig.  1  shown  by  Shettle  et  al.  (2010), 
however,  they  were  focused  on  inter-hemispheric  differences  and 
we  emphasize  interannual  variability.  In  addition,  in  contrast  to 
the  Shettle  et  al.  study,  we  have  removed  any  tidal  variability 
from  the  model-data  interpretation  by  sampling  NOGAPS-ALPHA 
at  the  locations  and  local  times  of  the  observations. 

The  top  panel  shows  the  NOGAPS  analyzed  temperatures  for 
all  the  SHIMMER  observations  in  July  for  the  50-58°N  latitude 
band.  The  bottom  panel  shows  a  histogram  of  temperatures  for 
those  SHIMMER  observations,  which  recorded  a  PMC.  Several 
points  are  immediately  evident.  First,  the  number  of  occurrences 
in  the  bottom  panel  is  much  less  than  in  the  top.  The  ratio  of  the 
area  under  the  bottom  curves  relative  to  the  top  curves  defines 
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Fig.  7.  Mean  temperatures  from  NOGAPS-ALPHA  for  the  month  of  July  corre¬ 
sponding  to  the  times  and  locations  of  the  SHIMMER  observations  in  the  50-58°N 
latitude  band. 
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Fig.  8.  (a)  Histogram  of  NOGAPS-ALPHA  analyzed  temperatures  at  all  the 
locations  sampled  by  SHIMMER  for  the  3  years  indicated,  (b)  Histogram  of 
NOGAPS-ALPHA  analyzed  temperatures  for  those  SHIMMER  observations,  which 
recorded  a  PMC. 


the  frequency  of  occurrence  for  SHIMMER  that  is  presented  in 
Table  1.  Second,  in  the  top  panel,  the  distribution  of  temperatures 
for  2007  is  clearly  shifted  to  greater  values  than  for  2008  and 
2009.  The  mean  temperature  for  the  2007  distribution  is  about 
160  K;  the  mean  for  2008  and  2009  is  155  and  154  K,  respectively. 
These  higher  temperatures  are  reflected  in  the  low  PMC  occur¬ 
rences  for  2007  in  Fig.  8b.  Third,  the  temperature  distributions  in 
the  bottom  panel,  those  associated  with  the  cloud  detections,  are 
clearly  to  the  left  (i.e.  colder)  than  the  overall  distributions  in  the 
top.  This  is  not  surprising.  What  is  somewhat  surprising,  although 
not  unprecedented,  are  two  features  of  the  bottom  panel.  First, 
the  temperatures  for  cloud  detections  for  the  three  years  differ. 
In  other  words,  there  is  not  one  representative  temperature  for  a 
PMC  from  one  year  to  the  next.  The  second  surprising  aspect 
of  Fig.  8b  is  that  there  are  considerable  cloud  detections  for 
temperatures  normally  thought  to  be  too  high  for  PMCs  (i.e. 
greater  than  150-155  I<).  This  relationship  is  also  particularly 
surprising  since  the  NOGAPS-ALPHA  analyzed  temperatures  and 
saturation  were  previously  shown  to  be  a  good  qualitative 
indicator  of  the  occurrence  of  bright  clouds  seen  by  SOFIE  in  the 
troughs  of  the  5-day  wave  in  late  summer  at  higher  latitudes 
(Nielsen  et  al.,  2010). 

This  wide  range  of  temperatures  associated  with  PMC  occur¬ 
rence  was  also  seen  by  Shettle  et  al.  (2010).  They  suggested  that 
variations  in  gravity  wave  activity  or  availability  of  nucleation 
sites  may  play  a  role  in  the  interpretation  of  PMC  occurrences 
between  hemispheres,  which  may  also  apply  to  our  interannual 
study.  In  addition,  there  are  several  other  points  worth  noting. 
First,  PMC  detections  for  limb-viewing  temperatures  greater  than 
the  presumed  frost  point  of  ~147I<  are  not  limited  to  MLS 
temperatures.  Stevens  et  al.  (2001)  reported  similar  observations 
using  CRISTA  data.  Since  CRISTA  temperatures  are  likely  more 
accurate  in  the  mesopause  region  than  MLS,  the  problem  may  not 
simply  be  due  to  MLS  bias.  Second,  it  is  clear  that  compared  with 
Shettle  et  al.,  there  are  many  more  SHIMMER  PMC  observations  at 
high  temperatures,  i.e.  160  K.  But  this  should  not  be  surprising 
since  as  a  limb  viewer,  SHIMMER  is  more  sensitive  to  the  dim 
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clouds  than  SBUV.  Thus  it  appears  that  the  few  clouds,  which 
Shettle  et  al.  (2010)  reported  for  T  >  150  K  were  only  a  fraction  of 
what  were  ultimately  present. 


4.  Microphysical  modeling 

Here  we  quantitatively  evaluate  the  implications  of  the  2007 
temperature  enhancement  and  its  implication  for  clouds.  Our 
approach  is  similar  to  that  presented  by  Stevens  et  al.  (2010).  This 
requires  producing  a  mean  temperature  field,  which  preserves  the 
local  time  information  so  that  a  comparison  can  be  made  at  the 
proper  time  of  day.  We  thus  averaged  together  the  124  6  hourly 
synoptic  samples  of  temperature,  H20,  horizontal  and  vertical 
winds  for  both  July  2007  and  2008  and  placed  them  on  a  local 
time  grid.  As  discussed  by  Stevens  et  al.  (2010)  this  removes  all 
longitudinal  variability  but  isolates  the  effects  of  the  migrating 
tides.  One  difference  from  Stevens  et  al.  (2010)  is  that  we  are  only 
interested  in  modeling  PMC  at  the  single  local  time  of  the  SOFIE 
measurement  (~23  h).  We  therefore  only  use  one  parcel  calcula¬ 
tion,  rather  than  the  24  used  previously  to  generate  a  complete 
diurnal  variation  at  hourly  intervals.  As  discussed  by  Stevens  et  al. 
(2010)  and  as  apparent  in  Figs.  7  and  8,  the  mean  temperatures  at 
SHIMMER  latitudes  are  too  high  to  properly  represent  the  condi¬ 
tions  under  which,  the  PMCs  occur  at  those  latitudes.  However, 
since  PMCs  are  a  near-constant  occurrence  at  SOFIE  latitudes,  the 
mean  temperatures  should  be  much  more  appropriate  and  we 
limit  our  analysis  here  to  the  SOFIE  conditions. 

As  discussed  in  Stevens  et  al.  (2010),  while  our  intent  is  to 
model  PMCs  coincident  with  SOFIE,  it  is  unrealistic  to  simply 
initialize  CARMA  from  a  “cold  start”  state  that  is  highly  super¬ 
saturated.  It  is  more  likely  that  the  clouds  grow  along  a  specific 
trajectory  according  to  the  prevailing  wind  field.  We  therefore 
follow  an  air  parcel  for  two  days  using  the  horizontal  winds  from 
NOGAPS-ALPHA  and  apply  a  constant  progressively  decreasing 


temperature  offset,  taken  from  Stevens  et  al.  (2010),  until  we 
reach  the  SOFIE  location  with  appropriate  temperature  and 
moisture  conditions.  Consistent  with  the  approach  of  Stevens 
et  al.,  we  use  one  water  vapor  profile  at  69  N  and  allow  the  water 
to  vertically  redistribute  through  sedimentation  and  sublimation. 
The  temperatures  and  vertical  winds  vary  along  the  path  as 
discussed  above  until  the  parcel  reaches  the  average  SOFIE 
latitude  in  July  of  68  N.  Other  inputs  such  as  the  rate  at  which 
the  model  cools  to  ambient  conditions,  the  eddy  diffusion,  the 
meteoric  smoke  distribution  and  the  expression  for  the  equili¬ 
brium  vapor  pressure  over  ice  are  the  same  as  used  by  Stevens 
et  al.  (2010). 

Fig.  9  shows  sample  output  from  the  trajectory  calculation  for 
both  2007  and  2008.  Fig.  9a  and  c  shows  the  temperature  along 
the  parcel  path  at  the  indicated  latitude  and  time  for  2007  and 
2008,  respectively.  Variability  of  these  temperatures,  for  example, 
the  warming  at  1 .8  days  in  Fig.  9c  is  most  likely  a  local  time  tidal 
effect  as  discussed  by  Stevens  et  al.  (2010).  Also  superimposed  on 
the  temperature  contours  in  Fig.  9a  and  c  is  a  solid  black  line, 
which  represents  the  altitude  of  the  calculated  cloud.  Focusing  on 
the  2008  case,  we  see  that  the  cloud  appears  to  form  at  85-86  km 
before  finally  settling  to  83  km  at  the  end  of  the  trajectory.  Berger 
and  von  Zahn  (2007)  have  addressed  the  question  of  at  what 
altitude  nucleation  begins.  They  suggested  that  cloud  formation 
begins  no  more  than  3  km  higher  than  the  final  observation  at 
69°N.  Our  results  are  generally  consistent  with  this  and  further 
suggest  that  it  is  the  temperature  increase  at  85-86  km,  not  at  the 
cloud  altitude,  which  is  responsible  for  the  interannual  cloud 
difference. 

The  growth  of  the  cloud  for  2008  is  indicated  in  Fig.  9d  with  a 
final  column  ice  mass  of  64  pg  m-2.  This  compares  well  with  the 
SOFIE  observation  of  69  pg  m-2  for  July  2008.  However,  the  2007 
calculation  yields  insignificant  clouds  (cf.  Fig.  9b).  This  is  consis¬ 
tent  with  our  suggestion  earlier  that  the  2007  temperatures  were 
marginal  for  cloud  growth  at  SOFIE  latitudes,  although  clearly 
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Fig.  9.  (a)  Evolution  of  temperature  along  the  2  day  parcel  trajectory  used  by  CARMA  for  the  2007  simulation,  (b)  associated  calculated  IWC,  (c)  same  as  (a)  but  for  2008 
and  (d)  same  as  (b)  but  for  2008.  The  lower  x-axis  shows  the  absolute  integrated  time  of  the  model  calculation;  the  upper  x-axis  shows  the  latitude  of  the  parcel.  The  black 
line  in  (a)  and  (c)  marks  the  altitude  of  the  calculated  ice  layer. 
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SOFIE  did  report  a  much  larger  IWC  in  2007  than  our  baseline 
calculation  indicates.  To  explore  this  further,  we  performed 
sensitivity  studies  where  we  arbitrary  assigned  various  constant 
temperature  biases  to  the  NOGAPS-ALPHA  temperature  field  that 
is  used  in  CARMA.  The  biases  tested  are  3  I<  or  less  and  within  the 
reported  systematic  uncertainties  of  the  MLS  and  SABER  tem¬ 
perature  retrievals  (Remsberg  et  al,  2008;  Schwartz  et  al,  2008) 
that  are  assimilated  by  NOGAPS-ALPHA.  The  results  are  shown 
in  Fig.  10. 

Fig.  10  shows  the  results  of  10  calculations  of  IWC,  5  for  each 
year.  Each  model  run  assumes  the  same  trajectory  and  cloud 
lifetime  (i.e.  the  time  over  which  the  cloud  is  assume  to  grow). 
However,  to  illustrate  the  sensitivity  of  our  results  to  possible 
uncertainties,  we  varied  the  temperature  in  each  case  by  a 
constant  offset  ranging  from  -3  to  +1  K.  Also  shown  in  the 
figure  are  the  IWC  observed  by  SOFIE  for  2007  and  2008  as 
the  horizontal  lines.  As  with  Fig.  9,  it  is  apparent  that  in  2007,  the 
conditions  for  ice  growth  were  marginal  at  best.  Using  the 
NOGAPS-ALPHA  temperatures  with  zero  assumed  temperature 
bias  yields  no  clouds.  However,  the  model  yields  increase  drama¬ 
tically  as  a  temperature  offset  is  introduced  from  -1  to  -3  K. 
This  may  reflect  the  fact  that  our  monthly  mean  atmospheric 
fields  are  averaging  in  periods  when  the  atmosphere  was  too 
warm,  thus  skewing  our  results.  The  2008  simulation  is  less 
sensitive  to  the  assumed  temperature  bias.  Using  NOGAPS-ALPHA 
with  no  bias  yields  a  calculated  IWC  that  is  within  ~20%  of  the 
SOFIE  observations  in  2008.  Assuming  progressively  colder  tem¬ 
peratures  increases  the  calculated  ice  by  a  factor  of  4.  The  figure 
clearly  shows  that  regardless  of  the  assumed  temperature  bias, 
the  model  predicts  greater  IWC  in  2008  than  in  2007,  which  is  in 
agreement  with  the  observations. 

In  addition  to  the  arguments  presented  by  Shettle  et  al.  (2010) 
concerning  possible  temperature  biases,  we  can  suggest  two 
other  reasons  why  a  negative  model  temperature  offset  might 
be  appropriate.  First,  in  2007,  not  all  the  SOFIE  observations 
recorded  PMC.  Only  92%  recorded  measurable  IWC.  The  mean 
temperature  for  that  92%  of  the  SOFIE  data  is  colder,  by  about 
0.5  K,  than  a  straight  average  over  all  the  data  and  our  NOGAPS- 
ALPHA  average  does  not  account  for  that.  In  2008,  this  is  not 
applicable  since  SOFIE  essentially  recorded  a  cloud  in  almost  100% 
of  the  occultations.  Second,  and  probably  more  significant,  is  that 
the  southward  yaw  of  SABER  occurs  in  mid-July.  We  have  found 
that  the  sudden  denial  of  SABER  to  the  NOGAPS-ALPHA  assimila¬ 
tion  seems  to  result  in  a  small  (1-4 1<)  jump  in  the  analyzed 
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Fig.  10.  Calculated  IWC  for  several  CARMA  model  simulations.  Each  simulation 
presents  the  end  result  of  a  2  day  calculation  as  shown  in  Fig.  8.  The  differences 
between  the  simulations  represent  various  constant  temperature  offsets  applied 
to  the  input  temperatures.  Results  for  2007  and  2008  conditions  are  presented. 
Also  shown  for  comparison  are  the  observed  IWC  measured  by  SOFIE  (as  the 
horizontal  lines). 


temperature  poleward  of  50°N  and  above  80  km  (not  shown).  By 
averaging  over  the  entire  month,  we  may  be  averaging  over 
temperature  fields,  which  are  slightly  offset.  Whether  this  is  an 
artifact  of  the  how  the  SABER  and  MLS  instrument  temperature 
biases  are  handled  in  the  assimilation  is  unclear;  however,  it 
serves  as  a  useful  reminder  that  the  analysis  depends  both  upon 
data,  which  have  uncertainties  and  upon  assumptions  about  how 
to  assimilate  the  data,  which  can  also  introduce  small  biases  that 
can  convey  more  significantly  to  the  simulated  microphysics. 


5.  Summary  and  conclusions 

The  present  work  has  provided  insights  related  to  interannual 
variability  in  PMCs  through  teleconnections.  First,  we  have 
explained  the  dominant  variation  in  NH  PMC  over  the  years 
2007-2009.  The  most  dramatic  variation  was  the  decrease  in  July 
2007  relative  to  2008  and  2009.  This  decrease  was  manifested  in 
several  consistent  ways.  Thus  we  saw  an  8%  decrease  in  the 
frequency  of  occurrence  in  clouds  seen  in  solar  occultation  by 
SOFIE  near  68°N,  a  factor  of  two  decrease  in  the  ice  water  content 
observed  by  SOFIE  and  a  factor  of  5-6  in  the  number  of  clouds 
seen  by  SHIMMER  at  sub-arctic  latitudes.  The  2007  PMC  reduc¬ 
tions  are  seen  to  be  linked  to  SH  winter  weather  as  first 
postulated  by  Becker  et  al.  (2004)  and  discussed  by  Karlsson 
et  al.  (2007,  2009b).  This  is  consistent  with  earlier  reports  of  such 
inter-hemispheric  teleconnections;  our  results  are  the  first  to  use 
a  global  analysis  of  the  stratosphere  and  mesosphere  with  a  cloud 
formation  model  to  document  these  effects. 

It  should  be  noted  that  it  remains  somewhat  unclear  exactly 
how  this  inter-hemispheric  teleconnection  operates.  An  early 
investigation  of  the  response  of  the  middle  atmosphere  to  a 
stratospheric  warming  perturbation  found  the  effects  to  be 
generally  limited  to  the  winter  hemisphere  (Garcia,  1987).  One 
limitation  of  this  work  may  have  been  the  use  of  linear  Rayleigh 
damping  to  describe  small-scale  frictional  dissipation.  More 
recent  work  with  models  that  utilize  a  more  complete  parame¬ 
terization  of  small-scale  gravity  waves  (Karlsson  et  al.,  2009a; 
Kornich  and  Becker,  2010)  have  found  that  the  anomalous 
circulations  generated  by  wintertime  wave-mean  flow  interac¬ 
tions  can  generate  temperature  anomalies  in  the  tropics,  which 
by  changing  the  thermal  wind  balance,  can  be  mediated  to 
anomalies  at  high  summer  latitudes.  Recent  observations  of  a 
link  between  the  QBO  and  inter-hemispheric  coupling  (Espy  et  al., 
2011)  are  also  consistent  with  these  ideas.  The  present  work 
certainly  supports  the  robustness  of  the  winter-summer  linkage, 
but  a  global  theoretical  explanation  for  this  linkage  is  beyond  the 
scope  of  the  paper. 

As  a  corollary  to  understanding  the  PMC  variability,  we  have 
compared  the  analyzed  temperatures  from  NOGAPS-ALPHA  with 
SOFIE.  While  there  are  10  K  differences  in  the  temperature  profiles 
between  86-90  km,  for  the  3  years  in  question,  the  interannual 
changes  displayed  by  NOGAPS-ALPHA  and  SOFIE  agree  very  well. 
Both  NOGAPS  and  SOFIE  show  that  at  the  SOFIE  latitudes,  the 
effect  of  the  temperature  anomaly  in  2007  was  much  larger  at 
altitudes  above  the  layer  where  the  ice  content  peaks  (i.e.  83  km). 
This  agreement  is  satisfying  confirmation  that  the  NOGAPS-ALPHA 
analysis  diagnoses  highly  vertically  structured  changes  in  upper 
mesospheric  temperature  and  links  them  to  larger  changes  else¬ 
where  in  the  atmosphere. 

The  question  of  the  non-uniqueness  of  the  temperature-cloud 
relation,  documented  by  Stevens  et  al.  (2001)  and  Shettle  et  al. 
(2010)  and  confirmed  here,  remains  a  puzzle.  Even  though  our 
comparison  with  bright  SOFIE  clouds  showed  good  qualitative 
agreement  between  temperature  fluctuations  and  cloud  occur¬ 
rence  (Nielsen  et  al.,  2010),  it  may  be  that  a  limb  viewing  satellite 
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has  too  coarse  horizontal  resolution  to  capture  local  temperature 
fluctuations  for  quantitative  correlations.  The  recent  result  that 
the  actual  ice  temperatures  are  lower  than  the  air  temperature 
(Petelina  and  Zasetsky,  2009;  Hervig  and  Gordley,  2010)  also 
suggests  that  the  limb  viewed  temperature  is  really  a  mix  of  both 
cloudy  (colder)  and  clear  (warmer)  airmasses.  However,  even 
ground  based  data  (Hoffner  et  al,  2003)  have  recorded  NLCs  at 
temperatures  as  warm  as  157  K.  Clearly  this  is  a  problem  deser¬ 
ving  of  future  study. 
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